Mass and temperature limits for blackbody radiation 
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A spherically symmetric distribution of classical blackbody radiation is considered, at conditions 
in which gravitational self-interaction effects become not negligible. Static solutions to Einstein field 
equations are searched for, for each choice of the assumed central energy density. Spherical cavities 
at thermodynamic equilibrium, i.e. filled with blackbody radiation, are then studied, in particular 
for what concerns the relation among the mass M of the ball of radiation contained in them and 
their temperature at center and at the boundary. For these cavities it is shown, in particular, that: 
i) there is no absolute limit to M as well to their central and boundary temperatures; ii) when radius 
R is fixed, however, limits exist both for mass and for boundary energy density pB- M < K Ms{R) 
and pB < Q/B?, with K = 0.493 and Q — 0.02718, dimensionless, and Ms{R) the Schwarzschild 
mass for that radius. Some implications of the existence and the magnitude of these limits are 
considered. Finally the radial profiles for entropy for these systems are studied, in their dependence 
on the mass (or central temperature) of the ball of radiation. 

PACS numbers: 



I. INTRODUCTION 



> 

^\ , The original motivation for this work is to explore the relation between entropy and gravity in a classical 

context (i.e. not quantum). The goal is to see if, under strong gravitational interaction circumstances, some 
] properties peculiar to black hole entropy (such as for example its area-scaling law), apparent in a semiclassical 
approach, can also be retrieved for conventional entropy in a purely classical scenario, in which the fields too, 
in addition to the background, are classical. Having this in mind, systems consisting entirely of blackbody 
radiation are chosen for investigation, with the hope that for them, both entropy can be readily calculated and, 
also, strong gravitational interaction circumstances can be reached. During the analysis of such systems we 
have found for them the existence of peculiar mass and temperature limits for each assigned radius R. 
^ ' The idea that under strong gravitational interaction conditions conventional entropy can acquire properties 

Pf] similar to black hole entropy is not new. In particular in Q it has been shown that the entropy of a spherically 
symmetric self-gravitating system in equilibrium becomes area scaling as the radius approaches the Schwarzschild 
value; in this limit the total entropy of the body is accounted for only by the entropy of the more external 
spherical shells. In .2] the entropy of a spherically symmetric system in equilibrium consisting of a perfect 
fluid with constant density has been explicitly calculated. The entropy has been found to loose its extensivity 
property and to reside more and more in the external spherical shells as the gravitational interaction goes 
stronger. The limit R = 2M [M is the mass, or mass-energy, as spelled out in Q, known also as ADM mass) 
cannot be reached in this case, as equilibrium is possible only if i? > (9/4)M (see, for instance, Q)- One 
would like to explicitly calculate entropy for equilibrium systems with other equations of state, possibly going 
to explore also the limit R ~ 2M. But it turns out that for all viable examples the radius of the system 
becomes infinite. This is a problem because in this case it becomes unclear what is intended with the concepts 
of "more external shells" and "mass" of the system. The recourse to different equations of state for different 
radii can give realistic models with finite size, but can hide the relation between entropy and gravity behind 
the circumstantial peculiarities of the constructed model. 

As well known (see for example 0) and we will recall later, a gravitationally self- interacting static distribution 
of blackbody radiation is in general precisely one of these infinite size systems. From a thermodynamical point 
of view this system is particularly appealing due to its simplicity, as its state is locally determined by only 
one thermodynamical observable, for example its temperature T . It is then tempting to use such a system to 
investigate the connection between entropy and gravity and we can imagine to circumvent the problem of its 
being of infinite size by forcing it to be adiabatically confined inside a finite size cavity. By applying a sufficiently 
high compression one could then expect the system to go towards its Schwarzschild limit. In principle it becomes 
then possible to study the behaviour of a pure classical system, this cavity filled with blackbody radiation (in 
our study all what is needed of blackbody radiation are purely classical properties, as we will see), in a strong 
gravitational regime. At these conditions, by the way, some quantum behaviour is expected to come on for this 
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system, as, of course, pair production or other effects P|. We decide however to study the system at a purely 
classical level, and to postpone any quantum consideration (pair production) at the end, once the characteristics 
classically expected are known. All this is the content of the following section. 

II. FIELD EQUATIONS AND THEIR SOLUTIONS 

We consider a perfect fluid consisting of a photon gas at thermodynamic equilibrium (somehow along the 
lines of 0, in which this same system is appoached, even if with different context and motivations) or, without 
quantum terminology, consisting of blackbody radiation. Its equation of state is 

P = p/3, (1) 

where p and p are respectively the energy density and pressure in the rest frame of the fluid element. Its 
stress-energy tensor is 

Tal3 = (P + p)UaUp + pga(3, (2) 

where Ua is the fluid element 4-velocity and gap is the metric tensor. In the case of static equilibrium con- 
figurations (necessarily with spherical symmetry), from Q and ((SJ Einstein field equations in Schwarzschild 
coordinates read 



(rr) (3) 

167r 

where the derivative symbol means derivative with respect to radial coordinate (area radius) r and where the 
metric is 

ds^ = -/(r) dt^ + h{r) dr^ + dil^ , (4) 

with / and h two non-negative functions of r, to be determined. 

For static gravitational fields, Tolman relation Q prescribes that T-y/— goo be constant throughout the system, 
where T is the temperature in the rest frame of the fluid clement. This means for us that 
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T{r)^f(r)= const. (5) 
For blackbody radiation, Stefan-Boltzmann classical result reads (in conventional units) 

P = \^T\ (6) 

with c the speed of light in vacuum and a the Stefan-Boltzmann constant. Applying this equation locally 
throughout our system (equivalence principle), from equation Q we have thus 

(J 

P{r)f'^{r) = const ^ (7) 

with C some constant with the dimensions of energy density (length^^ in geometrized units). Note that while 
the value of C depends on the chosen unit for length, the quantities p/C and ry/C are dimensionless. 
Making use of O, the first and the second of equations Q become 
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(8) 
(9) 

while the third, as can be verified (and as noticed in [^), is always identically satisfied if equations (|8I9|I are 
satisfied, as a confirmation of Tolman result. 

Before we start to solve these coupled differential equations, let us study some general characteristics which 
the solutions f(r) and h{r) should satisfy. h{r) can be put in the form 

= 132^' (10) 

(see for example Q) with m(r) given by 

m{r) = / 4Trr^p{r)df = — / 72?^^^' (H) 
Jo ^ Jo J v) 

that is, m{r) is the ADM mass inside r. 

For each allowed equilibrium configuration, in the limit r 

m{r) ~ ^7rrV(0), (12) 
so that m{r)/r 0. This implies h{Q) = 1 for each solution with p(0) finite. As regards /, from iQ we have 



p(0)/'(0) = ^, (13) 

SO that finite p(0) implies a non zero /(O). Taking the derivative of (|10|l and using and we have that 
in the limit r ^ 0, h'{r) ~ f 7r/9(0)r and f'{r) ~ r/(0)f 7rp(0), so that /i'(0) = and /'(O) = 0. 
We want now to solve equations (|8I9|I . We can write equation (as in '^) as 



, ifif + rf) 



and, substituting this expression of h in equation we obtain 



Cr3/ + 3r/3 ^ ' 

so that equations p4ll5|l are equivalent to the starting equations H8I9|I . Using, in addition, the previous results 
about the behaviour of / near the origin, the problem is thus reduced to find the solutions / of the second order 
equation 1)15(1 with initial conditions 



/(o) = vS^; /'(o) = o. (16) 

For each choice of the pair C, p(0) we have a solution /(r). Note however from equation that physical 
configurations p = p(r) depend only on the ratio C//^ and not on C and / separately. This implies that for any 
chosen p(0), the (infinite number of) solutions / that we find when C is varied, correspond always to the same 
physical configuration p = p{r). To fix p(0) corresponds to uniquely fix the physical configuration, irrespective 
of the separate values of C and /. Without loss of generality we can then study the solutions to equation 
(|15|l for a single arbitrary choice of the value of C (7^ 0); let us put for simplicity C = 1. All possible static 
configurations of our system are then in 1-1 correspondence with all the solutions / of equation 
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with initial conditions 



/(0) = V^= /'(0)=0, (18) 

where p(0) is any positive number and with h given by 114|l . Due to JB)), each specification of p(0) is equivalent 
to a specification of the central temperature Tc ■ 

Some symmetry considerations can now help to understand how the set of solutions / (and h) is arranged. 
Consider the transformation 

r f = Xr (19) 
f^f = >^f- (20) 

We have 

/(f) = A/(r7A) 

or 

fir) = A/(r/A) 

and then 

fir) = fir/X) 

and 

fir) = jf'ir/X) 

so that, if fir) is a given solution, /(r) is also a solution, as can be easily verified going through equation H17(l . 
fir) is such that /(O) = A/(0). On the other hand we know that, starting from a given solution /i, we can 
obtain all solutions if we change the initial conditions, /i(0) /a(0) = A/i(0), with A spanning all positive real 
values. Thus, in terms of an assigned solution fiir), a function fx is a solution if and only if 



/A(r)=A/i(r/A) (21) 



with A > 0. For the solutions h we obtain 



hxir) = hiir/X), (22) 

where hxir) and /ii(r) corresponds to fxir) and fiir) through ifTH) . 

Summing up, once we know the behaviour of solutions / and h for a given particular initial value /(O) = ai 
of /, arbitrarily fixed (> 0), we can readily get the corresponding behaviours of the solutions for every other 
assigned initial value a2- Let us study for simplicity the case /(O) — 1, that is the solutions to (|17|l with initial 
conditions 

/(O) - 1; /'(O) = 0. (23) 

From numerical integration of equation 1)17(1 and from equation ((14|l we obtain the solutions / and h reported 
in Figures n and 13 The structure of the corresponding equilibrium configuration can be seen in Figures 13 and 
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21 that respectively report p and 2m/r, as a function of r. Our choice C — 1 impHes that p and r are given here 
respectively in units of C and C"^^/^. From Q and (|21|l the properties of p for every other solution can readily 
be inferred. Note that p{r) approaches asymptotically as r — > oo. This implies that for this configuration 
and, by scaling property (|21|l . for every other solution, the radius of the system is infinite. 

Let us consider in some detail the functions h and 2m /r (Figures [21 and 0J. h starts from 1 at r = 0, and 
2m /r from 0; they reach their maxima, after which they go towards their own limiting values (5|, 



h{oo) = 7/4 (24) 

and 

(^) =3/7, (25) 

shown in the Figures, with damped oscillations around them, in which both the amplitudes of the deviations 
from the limiting values as well their derivatives at any order reduce asymptotically to 0. Analytically both 
limiting values can easily be verified. In fact consider for example that from H1U|I we have 



^ = (26) 

r h 

In equation © we see then that, at extremal points, that is when r has values r* such that h'{r*) — (that is 
(2m(r)/r)^^^. = 0), we have 



6 



. 04 h 



0.03- 



. 02 - 



0.01 - 



10 



FIG. 3: p versus r. 



2m{r* 



(27) 



As extremal values tend to (2m/r)ao when r oo, from equation Q considered in this limit, using (j26|l we 
have 



1 /2m 



= 1 - 



2m 
r 



2m 

r 



= 1-2 



2to 
r 



and thus equation H25() . as well the result (|24l) . Here we made use of the fact that in the limit r ^ oo, 
{r^ / f^y = implies /' — f/r. Note that equations (|26|l and (|22|) determine 2m/r for every solution; we see 
that 2m /r scales as h. From scaling property 122|) and from equations H18|) and (|20|) . the asymptotic behaviour 
of h and 2m/r when r ^ oo implies that in the limit of very large central density p{0) or very large Tc, 
h(r) ~ const — 7/4 and 2m{r)/r ~ const = 3/7, for every r, except the region of r very near to 0, in which h 
goes from 1 to 7/4 and 2m/r goes from to 3/7, with damped oscillations. 

As regards /, from Q we see that f'{r) > 0, Vr and from H27fl we have that in the limit r ^ oo, /(r) ~ -^7/3 r 
(see 0)- Due to equation this implies 



7 r h 

lim /A(r) = lim fx{r) = lim /A(r) = lim A/(r/A) = lim A - - = W - r, [r fixed) 

A— *u p(0)^oo io^oo A— *0 A^O OA V o 

and then 

1 3 1 

hm p(r) = lim p{r) = lim „ = — 5- {r fixed) (28) 

that is when the central density (central temperature) is allowed to increase without limit, the density (tem- 
perature) at each fixed radius r does not increase unlimitedly but approaches a limiting value proportional to 
r-2 (r-i/2). 

Finally, from scaling properties I|21I22|I . from initial conditions p8|l and from the behaviour of /(r), h{r) 
and p{r) near r = 0, shown in Figures and 13 we have that in the limit of very low central density, f{r), 
h(r) and p{r) are constant in a ever increasing range of r, so that the flat-space limit for constant low densities 
is recovered. 

To obtain a system with finite size we can imagine to enclose our radiation inside a perfectly insulating 
spherical cavity with area radius R, with its center coinciding with the center of the radiation distribution 
and to assume that thermodynamical equilibrium is reached, that is the temperature of the walls is equal to 
the temperature of radiation at r = i?. This will correspond to a certain central energy density pc for the 
cavity. Due to the local uniqueness of solutions to equation H17(l with initial conditions H18|l . this is equivalent 
to take a certain (physical) solution to equation (|17|l (that with p{0) = pc), to cancel all radiation outside 
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FIG. 4: 2m, I r versus r. 



r — R and to put an insulating spherical wall (of negligible thickness and mass) at r = i? exactly at the 
temperature the radiation has there. As regards energy density radial profile p(r), this means that there is a 
1-1 correspondence between all admitted cavities filled with blackbody radiation (that is every admitted central 
temperatures or boundary temperatures as well radii) and all possible p{r) distributions in finite ranges [0,i?], 
where p{r) — l/(87r/^(r)) and / is solution to (|17|l with initial conditions H18|l . 

Let us assume now, as allowed by BirkhofF theorem 0,01 that the Schwarzschild coordinates be chosen in 
such a way that outside the cavity exactly the Schwarzschild form of the metric obtains 



ds^ = -(1 - 2M/r) dt^ + ^- — - dr^ + r^dVt^. (29) 

^ ' ' 1 - 2M/r ^ ' 

The parameter M here is precisely the mass-energy or ADM mass as given by (|ll|l with r — R, that is 



M = / 47rrV(r)dr = - / -jj^dr, (30) 
Jo ^ Jo J v) 

This mass parameter M is the mass probed by Kepler-like experiments with gravitating bodies far away from 
the system. In our context it thus represents meaningfully the mass of the ball of radiation contained in the 
cavity. As the metric outside has been explicitly chosen, the values of the metric tensor components gap at 
the boundary are now fixed. Going at this boudary from inside, continuity of metric tensor components at the 
boundary of the cavity requires that 



f{R) = 1 - 2M/R (31) 
h{R) = 1/{1-2M/R). (32) 

Equations 1)10(1 . I|ll|l and H30|l show that h fulfills condition 1(32(1 . Making use of the freedom we have to choice 
C and / separately, provided the ratio C//^ be held fixed, we can always find a value of C such that the 
corresponding / fulfills the boundary condition 131(1 . In other words, once a single coordinate system is explicitly 
chosen, the arbitrariness in the choice of the constant C in equations (jH))-® is lost. In the case of a cavity filled 
with radiation, with central density pc, condition ((31(1 requires C to be 



^_P{R) _ {1-2M/R)^ 

where with /i we intend here the solution / for the case C = 1 and p(0) — pc- Each choice of central density 
/o(0) fixes the mass M within a given radius R and then separately C and /. 
From the results obtained so far some consequences can be drawn. 

A first consequence regards mass. We can infer it simply from the behaviour of the functions 2m(r)/r for 
equilibrium configurations, as detailed above. The non-vanishing asymptotic limit ((25(1 implies that there is no 
limit on the mass M; we have however that for each fixed R 
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FIG. 5: 2M/R versus z 



p{0)R^ {R is fixed). 



2M{R) 
R 



<K^ 0.493, 



(34) 



with K dimcnsionlcss not depending on R. The value of K can be read directly in Fig. 0] as the value of 
the maximum of the function 2m/r there reported. This means that, whichever is the energy density or the 
temperature of the blackbody radiation at equilibrium inside a cavity of radius R, its mass can never exceed 
a fraction K of the Schwarzschild mass Ms — R/2 for that radius. Let consider a cavity at its mass limit 
Miimit, for concreteness the configuration described in Fig. ^with area radius R equal to the extremal point. 
Suppose to supply it, in some unspecified manner, with some positive amount of energy AE, while mantaining 
fixed its radius. What we should obtain at equilibrium, if reachable, is simply a cavity with new conditions 
at the center and at the boundary and, by conservation of energy, with a mass M — Mumit + AE; but we 
already know that for the given radius an equilibrium configuration with such a mass value does not exist. We 
have therefore that no additional mass-energy can be stored at equilibrium in that given region of space if, in 
the absorption process, it is put in the form of blackbody radiation. For masses lower than the limit (|34|l . at 
least one equilibrium configuration exists, but also cases with 2, 3, 4, ... equilibrium configurations are allowed 
with different central temperatures; this can be readily understood taking into account the properties of the 
function 2m /r described above. Fig. |S1 permits to see, for a cavity with fixed area radius R at equilibrium, the 
dependence of its mass on the central density or temperature. It reports 2M/i? as a function of z = p{Q)R? . 
The plot is done in terms of z instead of p in order to obtain for 2M/ R an universal shape, i.e. not depending 
on R. The expression for z indicates how to rescale the abscissa to obtain, from the given plot, the actual 
function 2M{p{Q))/R for each fixed value of R. The explicit form of the plotted function, expressed in terms of 
the solution ft./(o)=i for initial conditions H23I) (the solution h presented in Fig. 12), is 



as can be easily derived from equation H26(l and from the scaling property (|22|l . In the Figure the maximum 
is, obviously, {2M / R)„iax — 0.493 and the extremal point z, is z = 0.439; this values can be used to calculate 
various quantities for the extremal equilibrium configuration. If, for example, the radius is the solar radius, 
R = Rq, we can calculate the extremal values of mass Mmax and central density p{Q)extr as follows: Mmax = 
{2M/R)max = 1.16-10^ Mq, with Mq the solar mass, and p{0)extr = z/R^^ = 0.91 • 10"^^ cm'^ to which 

corresponds, through (0, a temperature T^xtr = 1-95 • lO^^K. Note that this temperature gives kT^xtr = 1-68 
MeV, with k the Boltzmann constant. This implies that in the core of this system many photons are well above 
threshold for pair production and our assumption of a fluid consisting of blackbody radiation alone is no longer 
viable. Our results reside on two characteristics of blackbody radiation systems, namely the equation of state (Q) 
(and many systems besides blackbody radiation obey it, for example ideal Fermi gases in the extreme relativistic 
limit) and the peculiar energy-temperature relation ©. Only systems for which both these properties hold true, 
are describable in the way we have shown and the eventual existence of conditions could be investigated, for 
which this is the case for a mixture of blackbody photons and electrons and positrons at equilibrium through 
the reaction 77 ^ e+e~. On the other hand, as T^xtr oc when the radius is of the order of '-^ 100 Rq or 
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larger the contribution of pair production starts to be highly depressed and thus in any case our assumptions 
and results in principle permit to describe faithfully the system, not only at low densities but even up to its 
mass limit. As the mass limit (|34ll is very high, systems of this size can store in this manner mass energies of 
orders of 10^ Mq while remaining at equilibrium. Quite interestingly we can have thus physically allowed and 
easily (classically) describable systems, which are extremely compact, with lO^M© masses within a size of few 
hundreds of light travel seconds (i.e. a fraction of AU). 

In Fig. 13 we see that once the extremal central energy density or temperature are overcome, equilibrium 
configurations are still possible, at lower values of M. We have thus larger central temperatures, which give rise 
however to a lower total mass for the system. This is by no way a problem for energy conservation, obviously. 
The total energy of the system, given by M , contains in fact both the component of proper energy p and the 
(negative) component of gravitational potential energy fsj . Evidently, when the extremal temperature is 
overcome, the system settles down into a new equilibrium configuration with a stronger gravitational energy 
component. 

A second consequence regards energy density or temperature. Note that, on the basis of the correspondence 
mentioned above between cavities and solutions p{r) to equations l|3Jl, we immediately have the following. While 
there is no limit on the central density pc (or central temperature) of a given cavity at equilibrium (irrespective 
of being its size fixed or not), from (|28|l on the contrary when the size is fixed the temperature at its boundary 
(that is the temperature of the blackbody radiation coming out through a small hole in the walls of the cavity) is 
limited, whichever the central temperature is. Fig. [^reports energy density at the boundary of the fixed-radius 
cavity p{R) times as a function of the quantity z defined above. Here again we have an universal plot, not 
depending on R. Its analytic expression in terms of the solution //(o)=i corresponding to initial conditions (|23|l 
is 



//(o)=i(V87rz) 

as can be deduced from equation (|7|) and from the scaling property (|21|l . This dimensionless function has 
absolute maximum Q = 0.02718 (the extremal point is z = 0.1021) and goes to the asymptotic value 3/(567r) 
(with damped oscillations around it) when r — > oo as can be seen from H28|) . For a cavity of size R we have thus 
the following limiting density on the boundary 



PiR) < ^ (37) 

The limit depends on the size R of the cavity. For temperatures at the boundary of the cavity higher than 
those corresponding to this limit, thermodynamical equilibrium is no longer possible. If we perform a small hole 
in the cavity, the radiation emerging from it cannot have an energy density (or corresponding temperature or 
pressure) larger than this limit. Note that p{R) never vanishes, whichever is the {y^ 0) value of central density. 
Due to this implies that for our cavities, self-gravitating equilibrium configurations cannot exist. 

We study, finally, how the radial profiles for entropy density in the cavities depend on central density (or 
temperature) and then on the mass of the radiation inside it. As already said, this study is performed to 
see if the results P| on the variation of entropy radial profiles with mass for spherically symmetric bodies with 
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constant density can be extended/confirmed for other spherically symmetric equilibrium configurations. Instead 
to refer to entropy density in the proper orthonormal frame tT(r), following ^ we present our results in terms 
of, call it, entropy radial density s{r) defined by 



Attt^ sir) dr, (38) 

^0 

where S is the total entropy of the cavity and i?, as before, its radius. We see that s{r) is the specific contribution 
to S per unit area per radial slice (r, r + dr). Evidently the relation between s and cr is s(r)dr — a{£)d£, where 
di is the proper line element, so that (from (0}) 



di 

s = a— = crVh. (39) 
dr 



From first law of thermodynamics and assuming extensivity to hold locally (from equivalence principle) for the 
pertinent thcrmodynamical quantities, the Gibbs-Duhem relation Q 



p = Ta - p + (40) 

obtains, where /i is the chemical potential and n the particle (photon, in this case) number density in the rest 
frame of the fluid. Taking into accont that for photons n — 0, from and JHJl we obtain 

s^^^Vh^p'/'Vh. (41) 
From this equation and from scaling properties (|21I22() . we have that s scales as 



(42) 



where si is the entropy radial density for the reference configuration, corresponding to the functions /i and hi 
in equations H21I22|I . Putting x = r/R, with R assigned, we have 



sx{Rx)^ ^si{Rx/\). (43) 

This law express the scaling (|42|) in terms of the normalized radius x. 

Let us choose as radius of the cavity the radius Rq for which 2m /r in Fig. 01 is maximum. In order to 
compare entropy radial profiles when the central density is changed, Fig. [T] plots radial profiles of the function 
f{0)^^^ p^/^^/h for various choices of /(O). The radial coordinate is expressed as a; = r/Ro and p is intended 
again in units of C. The curves are with increasing f(0) upwards, respectively /(O) = 2/3, 1, 10, 100. This 
means that the second curve from below corresponds to the critical configuration, (for which the central density 
p{0)extr is such that M is maximum) and the curves above and below it are for central densities respectively 
lower and higher than p{0)extr- The Figure shows that at increasing central density the entropy distribution is 
more and more weighted around the center; further inspection of equation (|43|l gives that in the limit p(0) — *■ oo 
the plotted function f{0)^^^p^/'^y/h can be approximated as ~ const ■ [A/ {xRq + A) J'^/^, where A — /(O). Let us 
put together this Figure with Fig. If we start from low densities (or temperatures) we see that at increasing 
M the contribution to total entropy is more and more dominated by the most internal layers in r, at variance 
with what obtained, for another system, in where at increasing M the dominating role is gradually played 
by external layers. We see also that when however we are at densities a bit larger than p{0)extr the weight of 
external layers tends to increase with M, confirming in this case what found for the system studied in j^]- This 
will be reversed again after the next extremum for 2M/Rq (a minimum) will be encountered and so on. We 
see therefore that in our circumstances the dependence on mass of radial entropy profile is different and more 
complex than in the case examined in 0. 
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FIG. 7: f(0)^^^ p^^'^y^ versus x = r/Ro for various choices of central temperature, see text. The curves go upwards at 
decreasing central temperature. The second curve from bottom corresponds to the limit mass configuration. 



Our results can be summarized as follows. Considering insulated spherical cavities thermodynamically at 
equilibrium, i.e. filled with blackbody radiation, at energy densities sufficiently high to make not negligible the 
gravitational self-interaction of radiation, we determine their mechanical equilibrium configurations. We find 
that equilibrium is possible without any limit on mass M (ADM mass) of the contained radiation, provided the 
area radius R of the cavity is arbitrarily large. At fixed i?, on the contrary, a (very large) mass limit is present, 
which, if expressed in terms of the Schwarzschild mass Ms for that radius, is M{R) < KMs{R) with K = 0.493 
not depending on R. This implies that no more than this amount of mass energy can be stored as blackbody 
radiation in a cavity with assigned area radius R. Furthermore we find that equilibrium is possible for any 
value of central energy density (or temperature), irrespective of being the size of the cavity fixed or not. The 
energy density ps — p{R) at the boundary of the cavity can take any non-zero value if the size R of the cavity is 
arbitrary; it is however bounded if R is fixed. We find p{R) < Q/R? with Q = 0.02718 not depending on R. This 
can be interpreted as a limit on the density of the radiation emerging from a small hole in a cavity of assigned 
radius R. The fact that pB is never (except in the trivial case for which p{Q) = 0) implies that self-gravitating 
equilibrium configurations are never allowed. In the limit of very high central density (or temperature) we have 
an asymptotic situation in which both the quantities 2M/R and p{R) R^ no longer depends on R, being the 
first fixed to the value 3/7 and the second to 3/(567r). 

These results are obtained considering the radiation contained in the cavities from a purely classical point 
of view, that is neglecting any quantum process such as, for example, pair production by photons. This has 
been purposely done to study the radial profiles of conventional entropy density for classical fields under strong 
gravitational interaction conditions. For central temperatures not very high this description can be in any case 
safely adequate. It turns out furthermore that when the considered radii are roughly ~ 100 Rq or larger, 
these systems can reach even their mass limits without any significant appearance of pair production processes. 
Cavities with 100 Rq size (a fraction of AU) can be at equilibrium (i.e. not collapsing) even with masses of 
orders of IO'^Mq. 
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